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Abstract 



To explain the main features of the metal-insulator transition (MIT) in 
2D we suggest a simple model taking into account strongly localized states 
in the band tail of 2D conductivity band with a specific emphasize of a role 
of doubly-occupied states (upper Hubbard band (UHB)). The metallic be- 
havior of resistance is explained as result of activation of localized electrons 
to conductance band leading to a suppression of non-linear screening of the 
disorder potential. The magnetoresistance (MR) in the critical region is 
related to depopulation of double occupied localized states also leading to 
partial suppression of the nonlinear screening. The most informative data 
are related to nearly activated temperature dependence of MR in strongly 
insulating limit (which can be in particular reached from the metallic state 
in high enough fields). According to our model this behavior originates 
due to a lowering of a position of chemical potential in the UHB due to 
Zeeman splitting. We compare the theoretical predictions to the existing 
experimental data and demonstrate that the model explains such features 
of the 2D MIT as scaling behavior in the critical region, saturation of MR 
and H/T scaling of MR in the insulating limit. The quantitative analysis of 
MR in strongly insulating limit based on the model suggested leads to the 
values of ^-factors being in good agreement with known values for localized 
states in corresponding materials. 
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I. INTRODUCTION 



The problem of apparent metal-insulator transition in 2D structures observed for dif- 
ferent 2D systems including in particular Si MOSFETs (see e.g. |J, [|]) and GaAs-AlGaAs 
heterostructures (see e.g. ||, J3J, ||); the extent review of the situation is given in @) is 
still far from complete understanding. At the present time the situation is disputed be- 
tween the two possible scenarios. According to the first one the unusual metallic behavior 
which apparently contradicts the scaling theory of localization || is a real manifestation 
of some new physical mechanism. According to the second one, the metallic conductivity 
is still expected to be suppressed at low enough temperature while the behavior observed 
until now is related to the mechanisms of more conventional nature like temperature- 
dependent disorder M or magnetic-field driven disorder ||. Despite of various mecha- 
nisms proposed, according to the conclusion of a review paper J7j "while each of these 
is capable of explaining one or another part of experimental observations, none of them 
provide a comprehensive picture". 

The theoretical approaches developed until now were mostly concentrated on the 
metallic state and discussed an effect of disorder and electron-electron interactions on 
this state while to the best of our knowledge practically no attempts were made to look 
at the problem starting from the insulating limit. The role of localized states was to some 
extent discussed in the papers ||, 0, flOU , [flT] where an effect of traps in oxide layer 



or interface region of Si MOSFETs was considered. However these localized states were 
actually external with respect to the states of 2D conduction band and only played a role 
of additional scattering centers while we mean the states responsible for the 2D transport 
deep in the insulating limit. Then, the traps - although by no doubt important for Si 
MOSFETs - can be hardly considered for other systems exhibiting the 2D MIT behavior. 

To our opinion, the clue to the problem of the 2D MIT can be related to the MR 
in strongly insulating limit (which can be in particular reached from initially metallic 
states in high enough fields). As it is known, the strong positive MR in parallel magnetic 
fields leading to a suppression of metallic state in the "critical" region stays to be one the 
puzzles of the 2D MIT. We especially emphasize that strong positive MR persists deep 
in the insulating state and is exhibited by the systems which were perfectly insulating 
at H = [0 . The specific feature of MR in this regime is the fact that its temperature 
dependence is close to the Arrenius law even if for H = the system has exhibited variable 
range hopping of Efros-Shklovskii type (which was observed for Si MOSFETS [12| and re- 
type GaAs heterostructures ||13|| ). We believe that the renormalization approaches starting 
from the metallic state (see e.g. ||14j|), possibly explaining MIT and suppression of metallic 
phase by magnetic field, fail to explain this MR deep in insulating phase. 

As it was shown earlier JI5), the activated behavior of MR in hopping regime is a 
signature of the hopping over the states of the UHB, that is over the doubly-occupied 
localized states. Because of the on-site spin correlations the Hubbard energy acquires 
in this case the Zeeman addition gfisH (g being the g— factor) and thus magnetic-field 
dependent addition is a universal one (at least for g^sH » T). Basing on these facts, 
we suggested that the UHB plays an important role in the problem of 2D MIT |1| . Note 



that a possible role of the multiply-occupied electron states in 2D MIT was also discussed 
0- 



in 
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In what follows we suggest a model which starts from the localized limit of the 2D MIT 
rather than from the metallic one and thus is concentrated mainly on the "insulating" 
and "critical" regimes. We believe that in the systems exhibiting the 2D MIT there 
exists a tail of strongly localized states below the bottom of 2D conductance band with 
localization lengths a smaller than predicted by scaling theory of localization || starting 
from the metallic regime; localization length a as a function of energy for these deep states 
is expected to have a critical divergency with energy increase with a cut-off imposed by 
STL at some threshold energy. The highest of these states according to the predictions 
by Kamimura |TS] are the double-occupied ones. Namely, the energy dependence of the 
Hubbard energy leads to formation of the peak of the UHB close to the bottom of the 
conductance band [18| . The metallic behavior of resistance in our model is related (like in 
the model |J) to temperature-dependent disorder originated in our case due to activation 
of electrons from the strongly localized states to the conductance band which leads to 
partial undressing of the disorder potential. The MR in the metallic regime is expected to 
be due to redistribution of electrons between different doubly- and single occupied states 
(driven by the Zeeman addition to the energy of double-occupied states) which also leads 
to partial suppression of screening of the disorder potential. (The similar concepts were 
suggested in Refs. [HI , || where depopulation of filled traps to delocalized states due to 



an increase of energy of the traps in magnetic field was considered). This simple model 
is shown to explain - at least qualitatively - most of the features of the observed 2D MIT 
(including scaling of resistance in critical region, H/T scaling of MR in the insulating limit 
and saturation of MR at large H). Basing on the fact that in the strongly insulating limit 
the values of ^-factor are not expected to be renormalized by Fermi-liquid effects, we have 
compared some of existing MR data with predictions of our model. A good agreement 
between g-factors extracted according to our predictions and their known values gives to 
our opinion a strong support to our model. 



II. FORMULATION OF THE MODEL 

We assume that the 2D conductance band has a band tail of strongly localized states 
originating due to localization of carriers in the potential relief imposed by the disorder 
potential. We believe that the localization length for these states behaves as 

a oc (e m - e)~ v (1) 

where e m is some energy which for 3D would correspond to mobility edge while v is an 
index of localization length. 

Actually STL predicts for 2D an absence of any mobility edge and allows only some 
threshold values of conductance separating regimes of weak and strong localization. The 
localization length is expected to be exponentially dependent on the value of conductance 
in the weakly localized regime and to decrease when the conductance tends to the critical 
value mentioned above. Thus one may expect that for a given degree of disorder there is 
some minimal value of the localization length a available in the STL (of the order of the 
mean free path). 

In our model we will assume that the localization lengths of the strongly localized 
states described by Eq.jl] is smaller than a. This assumption corresponds to a cut-off of 
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the divergency in r.h.s. of Eq.|l| at some 



while for higher energies the situation can be controlled by the STL. In other words, at 
e m we have a crossover from strongly localized states described by Eq.[l], to the states 
described by STL. While the states in the vicinity of e m are still strongly localized, the 
localization length exponentially increases with an increase of e — e m . 

In what follows we will also assume that the difference e m — e m is small with respect 
to characteristic energy scales and can be neglected. It means that the width of the band 
of states strongly localized in accordance to STL is much smaller than the width of the 
band of strongly localized states described by Eq.|l|. Note that the behavior of this sort is 
expected just for the "new" 2D structures exhibiting MIT since these structures exhibit 
high mobilities and, correspondingly, the values of a for them are expected to be large. 
Having these considerations in mind in what follows we will still use the term "mobility 
edge" for the energy e m , 

As it is known, the systems exhibiting 2D MIT are characterized by strong electron- 
electron interactions. In our model we will directly take into account the intrastate 
electron-electron interactions which were shown in [18j to be most important. Follow- 



ing the concepts by Kamimura fl8fl , we take into account the energy dependence of the 
Hubbard energy U : 



e 2 



U ~ — (2) 
Ka 

where a is energy-dependent localization length while k is the dielectric constant. While 
for 3D it is known to be strongly dependent on localization length (through the wave- 
vector dependence) there are convincing arguments to assume that for 2D the value of 
K independent or almost independent of a and close to its bulk value; these arguments 
were given in the paper [2(J where the exponent in the relation k oc 



a 2 < 0.2 However we still consider a possibility of a 2 ^ and thus we have 

U oc (e m - e) u (3) 

where u > v. 

As it is shown in the Appendix 1, the energy dependence of the Hubbard energy has 
significant consequences for the form of the upper Hubbard band. Namely, for strong 
electron-electron repulsion this band is peaked near the energy e m |TB[ and the peak is 
characterized by a steep decrease of DOS for e < e_ ;C where e_ iC is some threshold energy 
(see Fig. 1). When \i = e m , that is at the point of metal- insulator transition (in the sense 
of strong localization) all states below the mobility edge (possibly except of the states 
deep in the band tail) are doubly-occupied. This fact and a presence of a peaked DOS of 
D~ band below e m allow us to conclude that the metal-insulator transition takes place in 
the D~ band and its presence significantly affects the features of the transition. 

We would like to note that the D~ states in our case, strictly speaking, are different 
from "standard" impurity D~ states similar to negative hydrogen ion. In the latter case 
the two electrons are in the field of a single positive elementary charge. For the "tail states" 
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it is not the case since the electrons are considered to be localized in the potential wells 
created by charge density outside of the 2D sheet like the charges of the interface traps in 
the case of silicon MOSFETs or spatial charge fluctuation as in the case of GaAS-AlGaAs 
heterostructures. In this case the localization centers are by no means of the Coulomb 
nature. Thus the double occupation of some localized state does not mean a creation 
of negatively charged Coulomb center as well as single-occupied state is not similar to 
neutral center. (Nevertheless for the simplicity we will still use the notations D° and D~ 
for singly- and doubly-occupied states.) As it was mentioned above, the occupation of the 
localized states by single electrons or by electron pairs for the external "probe" particle 
is reduced to the effect of non-linear screening of the initial disorder potential. 

As for the rest of electron-electron interactions, we believe that the latter can be taken 
into account as a nonlinear screening of the initial disorder potential by electrons trapped 
to these states with an account of electron-electron correlations. As a result, the form of 
the band tail and the position of the mobility edge are expected to depend on the electron 
concentration n. For small variations of electron density Sn linearizing the dependence in 
question one obtains 

5e m = >y5n (4) 

We believe that the parameter 7 is scaled with the magnitude of the disorder potential 
subjected to the nonlinear screening and, correspondingly can be estimated on the base 
of estimates of efficiency of nonlinear screening as 

e 2 
mi, 1 / 2 

Thus both the chemical potential and the mobility edge are expected to be sensitive to n 
and one has 

5{e m - fj) = (>y + g- 1 )5n (5) 

where g t is some total density of states at e = e m . Thus the critical value of n = n c 
corresponds to fi(n) = e m (n). 

We assume that the initial ("bare") disorder energy is comparable to electron-electron 
interaction energy e 2 n l J 2 / k which ensures effective non-linear screening. This factor arises 
in "new" systems with larger mobilities and weaker disorder. In systems with larger 
degree of disorder the latter - at least in the insulating limit - is less sensitive to electron 
occupation. Another feature of these systems is a larger spatial scale of the disorder 
potential which makes even bare localization length to be large enough. 

While the occupation of the localized states corresponds to non-linear screening of the 
disorder potential, the occupation of the states above e m drives the system into regime of 
linear screening according to the considerations given in Thus the metal- insulator 
transition is accompanied by a transition from nonlinear screening to linear one; however 
even in the metallic state - at least in the " critical" region - one expects that short-range 
harmonics of the disorder potential are still screened by localized electrons. 

It is interesting to note that in the critical region when g t can be considered to be 
equal to the 2D band DOS = 2m/ 'nh 2 , the ratio of the terms 7 and g^ 1 in the r.h.s. 
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of Eq.[5] is of the order of the interaction parameter r s = 2 1 / 2 me 2 / nhpp {pp being the 
Fermi momentum) and thus the metal-insulator transition is dominated by the shift of 
the mobility edge. 



III. TEMPERATURE BEHAVIOR OF RESISTANCE 

In our model we relate the temperature behavior of resistance in metallic state, in 
accordance to 0, to temperature dependent disorder. Namely, we assume that the ion- 
ization of strongly localized states below e m leads to additional scattering for the mobile 
electrons qualitatively in the same way as ionization of occupied donor state create an 
additional charged scatterer. Indeed, the nonlinear screening of the disorder potential by 
the localized states is reduced when the number of localized electrons is reduced. We ap- 
preciate that quantitatively the contribution of the states with large localization lengths 
is not expected to be similar to contribution of donor state efficiently separated from the 
other localized states, but believe that the qualitative similarity still exists. However the 
efficiency of the scattering by ionized localized state is expected to depend on the local- 
ization length of the state and in general to increase with decrease of a because of the 
larger momentum transfer. 

Thus the total "classical" resistance can be written as 



P = Po + Pi (6) 



where 



Pi = Pi(|) 9 exp(-(^) (7) 

Here the exponential describes activation and thus according to our considerations T s ~ 
fi — e m while the T-dependence of preexponential (chosen as some power function) is 
related to dependence of scattering efficiency on a and thus on the energy of the state. 
As it is seen, this expression coincides with one given in |3j for the traps below the Fermi 
level. 

Note that in principle the "residual" resistivity po is also temperature dependent due 
to activation-induced increase of concentration of mobile carriers. To take this dependence 
into account we apply the Drude-like relations 

11 

Po = , Pi 



where l//io an d 1/pi are the inverse mobilities controlled by the "residual disor- 
der" and "temperature-dependent disorder", respectively. With a simplest assumption 
d(ni 1 )/dn = jl^ 1 = const (that is with a neglect of a dependence of scattering efficiency 
on a) one easily obtains 

pan-m = (- UN- 1 ) o» 

\n m o n m (l ) J \ H po/ 

where n m0 = n m (T = 0). 
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IV. MAGNETORESISTANCE 



The simplest situation corresponds to dielectric limit (/x < e m ) with an additional 
assumption go(fi) > <?- (/■*)• In this case (see Fig.l,b and also Appendix 2) the main 
effect of the magnetic field is related to the shift of \x in the D~ band with magnetic field 
increase. For T = this shift is equal to Sfi(H) = HsgH (where the value of g corresponds 
to localized states) while for finite temperatures the position of the chemical potential is 
described as (see ||T5| ) 

^f = Mcosh^) (10) 

In the regime considered the dominant contribution to conductivity is related to hops 
within the D~ band from doubly occupied states with e = fi since the hops with spin flips 
from the single occupied states to single occupied states are not allowed (as in the case 
of mechanism considered by Kamimura |JL8||). 

We will begin with the case when there is no Coulomb gap at the Fermi level that 
is when the interstate Coulomb interactions are neglected. Note that for the hopping in 
the band tail with small DOS the situation crucially depends on the relation between the 
width of effective hopping band Se 

0£ — ~, ; — ; m exp 7r - L - 11 

G? c (a/2) 2 )V3 3£ 2 ^ J 

estimated for a standard VRH way for the value of DOS equal to g — const = g(/i) (let us 
restrict ourselves by the VRH of the Mott type) and and the energy scale E\ characterizing 
energy dependence of DOS near the Fermi level. If be < S\ then one deals with a standard 
Mott-type hopping estimated for the corresponding density of states being approximately 
symmetric around the chemical potential. 

However when the opposite unequality holds which can be rewritten as 

^^>21og(fc £l ) 1/2 ^) (12) 

the hopping is expected to be dominated by electron hops from the Fermi level to empty 
states high above the Fermi level where the density of states drastically increases. 

We believe that it is legitimate to consider this sort of hopping in a standard for 
the VRH way, that is optimizing the hopping probability with respect to the typical 
hopping distance ~ N~ 1 / 2 (where N £ is a density of sites with energies less than e) and 
to the activation energy e — //. Here we consider a case of the Gaussian tail while the 
generalization for the exponential tail is straightforward. Taking into account that 



at Fa ' { £ c-ey 
N £ = de g c exp ^ ~ g c e x exp 



p 2 
-1 



(where the integral is considered to be dominated by the upper limit) we obtain the 
following condition for the activation energy e*\ 
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£l r, / £ l/_ . M/2 4 ^ , ln „ ei 112 



Since, at it is seen, S\/{e c — £*) — 1 we will neglect the last term in []. Correspondingly, 
the temperature behavior of conductivity is expected to obey the law 

Sc-^-^tlog^^) 1 /^)]^ 
er oc exp )- (14) 

As it is seen, the activation energy logarithmically depends on temperature and slightly 
decreases with temperature decrease. The difference with respect to a pure Arrenius law 
originates due to a finite density of states in the tail. 

If one includes the Coulomb gap into considerations, it would lead to a quadratic gap 
near the Fermi level, but the strong energy dependence of "bare" density of states does 
not allow for this gap to develop up to its nominal value according to the value of n(/j.) 
due to a cut-off at energies ~ S\ controlling the DOS decay. 

Now let us turn to the magnetic field dependence which is the most important for us. 
As it is clearly seen, for strong magnetic fields the position of the Fermi level is significantly 
lowered which leads to a dramatic decrease of DOS at Fermi level. In particular, the 
following scenario is possible. If at H = the DOS at the Fermi level is large enough 
(see Fig.l,a) for the unequality opposite to the one of Eq.[l2| to obey, the standard VRH 
of Mott type or of Efros-Shklovskii type is observed. In particular it is expected for the 
case when initial position of the Fermi level corresponds to the peak region of D~ band. 
However an application of magnetic field shift the position of D~ band with respect to 
the Fermi level and finally can lead to nearly activated behavior according to Eq.|14] when 
fx B gH > e x . 

One notes that this hopping via the states of D~ band actually competes with hopping 
over (spin-polarized) D° band, that is with hopping from occupied to empty states. Cor- 
respondingly, when magnetic field-dependent contribution of D~ band becomes smaller 
than the contribution of hopping over the D° states, MR saturates. 

If we start from the metallic regime (see Fig. 2, a) the variation of e m with magnetic 
field is of principal importance. Moreover, one expects the MR in this case to be mostly 
related to the variation in question. Note that such a behavior is in agreement with the 
ideas of "magnetic field driven disorder" formulated in ||. An important factor is related 
to the fact that the metallic state in principle corresponds to linear screening while the 
spin polarization of the mobile electrons is not expected to produce a strong effect on the 



disorder potential (see e.g. [23]. However we believe that at least in the critical region 
the linear screening still co-exist with the nonlinear one and, in particular, the localized 
D~ states well below e m with relatively small localization lengths a < n^, 1 ! 2 can not be 
efficiently screened by the mobile carriers. 

Thus the devastation of the D~ states corresponds to blocking of the nonlinear screen- 
ing which is similar to creation of new scattering centers. So the polarization of localized 
states can lead to significant shifts of e m (see Fig.2,b). The more detailed arguments are 
given in the Appendix 2. 

Estimating the variation of the electron concentration due to a suppression of the D~ 
states with magnetic field increase as 

5n ~ g_(fi B g5H) 
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and applying the arguments similar to those leading to Eq. [5] one obtains an estimate for 
the shift of the mobility edge as 



S[e m - aO — ig-^BQbH (15) 

So one concludes that an increase of magnetic field can lead to a change of a sign of the 
l.h.s. of EqfTS] from negative to positive which corresponds to a suppression of metallic 
state. 

One concludes that in this regime the saturation of the MR corresponds to stabilization 
of e m which corresponds to total polarization of the electron system including the localized 
states. The saturation field H sat corresponds to the lifting of the bottom of D~ band to 
chemical potential. 

As for the final state of the system for H > H sat , it obviously depends on the initial 
electron concentration. If it is large enough, in the final state the chemical potential is still 
larger than e m and the system persists in the metallic state. Moreover, we believe that 
deeply in the metallic state the linear screening dominates and the MR can be related 



to the mechanisms of a sort of considered in p3fl . Thus we do not intend to analyze this 
regime in detail. 



V. DISCUSSION 

Now let us compare predictions of our model with existing experimental data. We 
shall start with temperature behavior of resistance. According to our Eq.|9| the sign of 
temperature coefficient of resistance depends on the sign of second bracket in Eq.|9|. More- 
over, Eq.|9] describes "metal-insulator transition" within purely classical regime. Indeed, 
for 

n m o = fi/no = n mc 

p(T) = const. For smaller values of n m the temperature behavior of resistance is dom- 
inated by activation of electrons from localized states to delocalized ones (in accordance 



to observations reported in [] 1 



p oc exp( £m T /X ) (16) 

One notes that even for purely classical description our model exhibits with exponential 
accuracy scaling of resistance behavior p(T,n s ) which is symmetric with respect to "crit- 
ical" curve with e m — p = if one assumes a constant DOS at the "critical region" since 
in the latter case (see Appendix 2) 

n s -n c = g t (p - e m ) 

This behavior is in agreement with experiment (see e.g. J7|). For strongly localized regime 
when e m — p is large enough the system is in hopping regime. While deep in this regime 
variable range hopping is expected which is in agreement with experimental data (see e.g. 

E)- 
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Actually the quantum effects tend the system to insulator for T — > since the states 
above the cut-off energy e m discussed in the Section 2 are described by STL. The cor- 
responding scenario is considered in detail in |J and thus the metal-insulator transition 
discussed above is not a true one. 

Now let us discuss the MR. As we have mentioned above, we are mostly interested in 
the insulating limit because, to our opinion, it can give more information about the system 
properties than the metallic limit. Thus we would like first to compare our predictions 
concerning the MR in strongly insulating limit with existing experimental data. 

As it was reported (see e.g. |7j], citeKlapwijkl, ), the high-field insulator regime typi- 
cally exhibits nearly activated temperature dependence of the magnetoconductivity. We 
shall start from experiment reported in the paper |H| where a delta-doped GaAs/AlGaAs 
heterostructure was studied. Although strictly speaking this system did not exhibited 2D 
MIT, it demonstrated strong positive MR starting from the hopping regime. A specific 
feature was a gradual transition from the variable range Efros-Shklovskii hopping to nearly 
activated hopping. In our previous paper |16] we have tried to fit the experimental curves 
according to the pure Arrenius law ascribing the deviations from experimental behavior 
to a neglect of the non-zero DOS in the band tail. 

Now we are going to compare the experimental data of Ref. |I3J with our present 



theory taking into account a Gaussian tail of the D band. At the Fig. 3 we have 
plotted theoretical curves corresponding to our Eq.[L4] calculated for e c — fi(H = 0) = 
0.234meV(H = 8T), 0, 219mev(# = 6T), e x = O.UmeV , g a 2 = 0.93meV-\ g = 0.12. 
As it is seen, there is a good agreement between experimental data and our theoretical 
prediction which supports our present interpretation. 

The value of the effective g-factor is about 4 times lower than the handbook values for 
GaAs. However, one has in mind that we deal with AlGaAs-GaAs heterostructure rather 
than with the bulk GaAs. The ^-factor values for AlGaAs quantum wells were calculated 



theoretically |2"4 |. It was shown that due to the fact that in such structures one has a 
mixture of GaAs states (for which the ^-factor ~ —0.45 is negative) and AlGaAs states 
(where g -factor is positive) the effective ^-factor depends on the well width. In the case 
of gated heterostructure under discussion, the effective width of the potential well is in 
particular controlled by the gate voltage V g and thus the effective (/-factor is expected to 
depend on V g . The corresponding behavior was considered in some detail in our previous 



paper |16[ and was shown to be in agreement with the experimental data of |L3 



In addition, we have compared our results with MR data by Simonian et al. [Q. The 
data obtained for Si MOSFETs demonstrated a suppression of metallic state in strong par- 
allel magnetic fields and transition to strongly localized state. Since our model gives most 
definite quantitative predictions for insulating limit, we analyzed the data on temperature 
behavior of MR reported in || for the strongest magnetic fields (1.4 and 1.2 T). On Fig.4 
we present both experimental data and theoretical curves calculated on the base of Eq.|H| 
with the values of the parameters e c — fi = 0.147meV(H = 1.4T), 0.12QmeV(H = 1.2T) 
Ex = O.lmeV g^a 2 = \.2meV~ x , g = 1.7. As it is seen, there is a good agreement be- 
tween theory and experiment. The value of g factor estimated in this way is close to the 
textbook values for Si. Note that in the insulating limit Fermi-liquid renormalization of 
g-f actor is irrelevant. 

Following the same ideas we have also analyzed the MR in the insulating limit for the 
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p-type GaAs-AlGaAs heterostructures exhibiting the metal-insulator transition ||. The 
extracted value of g-factor has appeared to be equal to 0.025-0.04 which is in agreement 
with the fact that the expected value of g factor for heavy holes is zero and the effect is 
related to admixture of light hole states. 

The paper reported a suppression of the metallic phase at strong magnetic field 
followed by a saturation of resistance, the behavior exhibited H/T scaling while the 
most representative data corresponded to the dielectric side of MIT. The corresponding 
behavior is in agreement with our Eq.pl] predicting the scaling for the MR activation 
exponent. 

An important feature of the MR in the systems exhibiting 2 D MIT is its saturation at 
high fields. According to our model, the nature of this effect can be different for different 
regimes. 

For the systems initially being in the strongly insulating limit (when the disorder 
potential is not significantly affected by magnetic field), the saturation field H sat is related 
to a competition between contributions of D~ and D° subbands and thus depends mostly 
on the forms of the subbands. 

For the systems initially being in the "critical" region (of no matter on what side of 
the transition) the initial evolution of the MR is related to shift of e m , although at high 
fields a competition between hopping processes in D~ and D° bands is still possible. In 
both regimes mentioned above the value of H sat is expected to be sample-dependent. Note 



that such a behavior was reported in [I I 



Although from the very beginning we have not intended to discuss properties of the 
systems where initial electron density was significantly larger than the critical one, we 
would like to give some comments for this regime as well. In particular, even for the 
systems being initially deep in the metallic regime but still containing a significant density 
of the localized states the value H sat is not expected to be universal since the form of 
D~ band is obviously sample-dependent. Moreover, the possible exchange of electrons 
between the conducting band and the states spatially separated from the conducting 
paths (discussed in Appendix 2) does not allow to relate the total electron concentration 
to position of the chemical potential in conducting band. In addition, one also notes that 
the value of g-factor for mobile electrons is renormalized due to exchange interactions 
and is concentration-dependent; in any case its value is different from the corresponding 
value for localized states. Thus we believe that the saturation field itself for the metallic 
regime does not give a proper information concerning to say the value of ^-factor. This 
conclusion is especially important for the samples with concentrations close to the critical 
one. 

Note that if H sat corresponds to the crossover from "metallic" to "insulating" state 
then - at least for magnetic fields close to H sat - the chemical potential crosses the peak 
of the D~ band. This factor can explain the enhanced DOS at the Fermi level for the 
"critical" state (with respect to the standard one for the ideal 2D conduction band) 
observed in [ETJ. 

We would like also to note that the "final" insulating state achieved in strong magnetic 
fields is different from the insulating state existing for H = at small enough electron 
concentrations; in particular, the latter situation is characterized by a presence of a peak 
of D~ band close to e m which is absent for the former one. Correspondingly, the electrons 
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in metallic state close to the "critical" state for H = are expected to suffer additional 
mechanism of inelastic scattering related to activation of D~ states which is absent for 
spin-polarized situation. This latter factor can enhance a role of weak localization for 
H > H sat . In any case, according to convincing calculations of f22|, the behavior of the 



"critical" states is strongly sensitive to the system parameters which - with an account of 
the factor mentioned above - explains a difference in behavior of the " critical states" for 
the cases of H = and H > H sat reported in [FZ1[] . 

We believe that the surprising evolution of temperature coefficient of resistance in 
course of magnetic field increase (" metallic-like" in H = and H > H sat and " insulating- 
like" for some intermediate fields close to H sat ) reported in |12[ can be related to evolution 
of a position of D~ band with respect to the Fermi level. Indeed, one can expect that the 
presence of quasi-localized states at the Fermi level at the moment when the peak of the 
D~ band coincides with /i can significantly affect both classical and quantum contribution 
to resistance. 



VI. SUMMARY 



We have demonstrated that the main features of the 2D MIT apparently observed 
recently in different systems can be explained within the framework of the simple two- 
band model. The latter is based on the assumption that the transition-like behavior takes 
place in the disorder- induced "tail" of the 2D conductance band where the lower, strongly- 
localized states are not described by scaling theory of localization starting from the large 
conductance limit. An important role among these states is played by doubly-occupied 
states coexisting in general with the single-occupied ones. The role of the D~ states is 
emphasized by the fact that as a result of energy dependence of the Hubbard energy their 
DOS is peaked near the bottom of the 2D conductance band. We also emphasize a role of 
non- linear screening of the disorder potential by localized electrons. The "metallic-like" 
temperature behavior of conductivity is related to activation of the localized carriers to 
the conductance band leading to undressing of the disorder potential due to the partial 
suppression of the non-linear screening. The strong positive MR on the dielectric side is 
related to a suppression of activated hopping to the peak of D~ band as a result of lowering 
of chemical potential with respect to D~ band. The suppression of the metallic state in 
strong magnetic field is explained as a result of a suppression of doubly-occupied localized 
states participating in the non-linear screening of the disorder potential. The comparison 
of the theoretical predictions with existing experimental data exhibits at least qualitative 
agreement while the values of ^-factors extracted from the MR data for strongly insulating 
limit according to our model are in quantitative agreement with the values known for the 
localized states in corresponding materials. 
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VIII. APPENDIX 1. FORM OF THE D- BAND 



The energy of the doubly-occupied state e_ is obviously given as 



e + U ( £ -^^r (17) 

Em 



Correspondingly, the DOS in the D band is given as [fL8 



\de 



where e is the energy of single-occupied state while the functional dependence of £q(e- 
is given by a solution of Eq.|T7|. 

As it is seen, g_ diverges at some Eq^ corresponding to vanishing of the derivative in 
Eq.|TJ: 

e m = uU ( £ -^^r- 1 (19) 

This is related to the fact that the function Eq(e-) is a double- valued one and vanishing 
of the derivative in question corresponds to a coincidence of these solutions. One branch 
of the solutions corresponds to the single-occupied states above e ,c while another one - 
to the states below £ 0>c . 

Then, one concludes that if 

^■m Eq c — E m I I 
\ Uq / 

is small enough with respect to the depth of the band tail, the lowest localized states 
allow only single-electron occupation because the Hubbard energy for these states appear 
to be too high. 

The critical value of the energy of double-occupied state e_ jC is related to £o,c by Eq jl7|. 
To estimate the energy dependence of g_ near the critical value of E- we take into account 
that for £ close to £o,c 

8e- d 2 E- , , . 

_._( £o - £o ,J (20) 

To rewrite this equation in terms of e_ we take into account that 

e -- £ -'-W eo=ee + 2^eJ {£ °- £ ^ (21) 
Having in mind that the first derivative vanishes, we are left with 

/-■ o2 \ -1/2 

(eo-e , c )^ l~j (£--£-,c) 1/2 (22) 
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where 

d 2 e. u- 1 fuU \ ll{u - l) 



9Eq £ m 

Correspondingly, we have for e_ close to e_ iC : 

_ 2 1 / 2 g>(£ (£_)) £ m ^ 1/2/ , £ m q/ 2 ( m -i) /„„s 
[£- — £-,c) 1 U—l uUq 

Note that in realistic situation the sharp edge of the DOS singularity is expected to 
be smeared in some way. In particular, the local fluctuations of the parameter Uq are 
expected due to local fluctuations of the dielectric function. It is natural to assume that 
the fluctuations of Uq correspond to Gaussian distribution with some variance 5Uq. As it 
can be easily followed, these fluctuations give rise to Gaussian fluctuations of the energy 
e_ jC with the variance 

e x ~ SU — 

U — 1 V Uq J 

Averaging the form-factor function of Eq.|23| 

1 



(e_ - 



\l/2 



with respect to the fluctuations in question one easily obtains that the averaged form- 
factor has a form 

((^-£_, c ) + fe_r 1/2 (e_>e_, c 
( £iexp _(£_Z£z£)2) £ _< £ _ (C (24) 

As a result one expects the D~ band to have its own tail corresponding to non-zero but 
quickly decreasing DOS at E- < £_ >c : 

g_~g c exp-(- — -f (25) 

where we have omitted indices — at energy terms implying that the function g-(s) has 
energy argument corresponding to D~ band. 



IX. APPENDIX 2. OCCUPATION OF D~ BAND 

Now let us discuss the details of the occupation of D~ band with an increase of the 
concentration. If one neglects the states in the tail of D~ band, the occupation of this 
band starts when the chemical potential touches its bottom, i.e. when \i = e_ jC . The 
double-occupied states first appear in the band of single-occupied states near the energy 
£o,c symmetrically with respect to this value. Further increase of concentration leads to 
rise of the position of chemical potential in the D~ band and to a broadening of the band 
of double occupied states in D° band. Finally \i reaches the value e m and for /x — > e m 
the gap between empty states and doubly-occupied states in terms of energies Eq tends to 
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zero. In other words, at the moment of metal-insulator transition (in the sense of strong 
localization) all states below the mobility edge are doubly-occupied (possibly except of 
the states situated deep enough). Actually one should have in mind that according the 
considerations given above an increase of concentration affects not only a position of ^, 
but a position of e m as well and thus the picture is more complex one than a simple 
occupation of a band with a given density of states). 

Formally at low temperatures one has for the chemical potential the following equation: 

n = J"de{g + g-.), (26) 

where g correspond to DOS of single-occupied states in the band tail, g_ - to DOS of 
doubly occupied states discussed above. Since according to our considerations mentioned 
above the very form of DOS in the band tail in principle depends on n, the latter being 
given by a solution of the general Hubbard Hamiltonian including both intrastate and 
interstate interactions, the problem of calculation of \x for a given n is extremely difficult. 

However for the small variations of n near the critical value n c one expects the lin- 
earization procedure of the sort of Eq.^| is possible. 

Another problem of extreme importance for us is related to variation of \i with a 
variation of magnetic field. The presence of the field leads to the Zeeman addition HbqH 



in the equation [L7|. Since this addition is not energy-dependent, it does not change a value 
of £o,c discussed above, however shifts the bottom of D~ band £_ jC towards the mobility 
edge corresponding to majority spins. Correspondingly, the number of doubly-occupied 
states in D° band around the energy e ,c decreases. Finally the value e_ jC coincides with 
e m which obviously corresponds to 

e m - £_, c (# = 0) = n B gH (27) 

After that moment only single-occupied states exist below e m . Recall that the value of 
e m is considered to be related to majority spins while for minorities it is enhanced by the 
Zeeman splitting. Thus for the magnetic field given by Eq.^7] near the mobility edge for 
majority spins the mobile electrons with majority spins co-exist at the same energy with 
localized D~ states which in particular gives additional mechanism of resonant exchange 
scattering. 

As for the position of the Fermi level, for a fixed number of electrons it should be 
calculated with an account of a balance between the number of electrons in single- and 
double occupied states. For the system deep in the insulating state such analysis was 
made in []15| with an assumption that the DOS of the double occupied states at the Fermi 
level is much less than DOS for single-occupied states. In this case the position of \x in 
the D° band is fixed while its position in D~ band is lowered according to considerations 
given above. 

The picture is not as simple if we start from metallic state. If we would deal only with 
the band tail states, depopulation of double occupied states (including the localized D~ 
states) would inevitably lead to rising of the chemical potential for the majority electrons. 
Note that for conducting (delocalized) electrons the value of g- factor can be renormalized 
due to Fermi-liquid effects (see e.g. P5| . ) 

Thus these simple considerations concerning the evolution of chemical potential with 
increase of magnetic field can not explain the strong positive MR in the metallic state since 
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the rising of chemical potential for spin-polarized electrons is not expected to decrease 
the conductivity significantly. To our opinion, the MR can be related to evolution of the 
position of the "mobility edge" e m with increase of field. This argument is supported by 
the estimate given in Section 2 which shows that - at least at the critical region - the 
evolution of e m dominates the evolution of /J by a factor r s » 1. 

Two possible scenarios can be discussed. First, even if one restricts himself by the 
band tail states he can expect that the devastation of the doubly-occupied D~ states 
(that is redistribution of electrons between double - and single occupied localized states) 
lead to enhancement of the disorder potential and to increase of e m . Indeed, one has in 
mind that occupation of the localized states by electrons corresponds to screening of this 
potential and the screening is expected to be more effective if there exists an additional 
choice between double and single occupation (lowering the total energy of the system). 
According to the arguments given above the evolution of e m dominates over the inevitable 
increase of /x related to the devastation of D~ states. 

Another possibility corresponds to a presence of the subband of strongly localized 
states assumed to coexist with the 2D band states (including the band tail states), but 
being spatially separated from the band states. These localized states are expected to have 
a mobility edge much larger higher than e m and much larger Hubbard energies preventing 
the double occupation of these states. Such an assumption is a realistic one since for 
small electron concentrations the percolation character of electron transport is expected 
due to spatial variation of e m due to large-scale potential In this case the electrons 



from D~ states are in particular redistributed to these states decreasing the density of 
electrons in the band tail. Thus the nonlinear screening of localized states becomes more 
effective while the rising of chemical potential is suppressed. 

As a result, the difference n(H) — e m (H) can at some H = H c change its sign which 
corresponds to a crossover to insulating limit. 

Although a possible shift in e m with increase of carrier concentration was disputed 
in [[n| we would like to note that the authors based their considerations on an absence 
of significant shifts of e m for small variations of electron concentration around its critical 
value. To our opinion, this argument does not hold for our picture since the redistribution 
of electrons due to spin polarization is expected to be pronounced in contrast to the case 
of small variation of n,. 
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Figure captions 

Fig.l,a The structure of the subbands of double-occupied (left) and single-occupied 
(right) localized states for the electron occupation corresponding to insulating limit in the 
absence of magnetic field. Hopping conductivity is controlled by the states of D~ band. 

Fig.l,b The occupation of the subbands depicted in Fig.l,a for high magnetic field. 
Hopping is dominated by activation to the peak of D~ band (e a = e c — fi) 

Fig. 2, a The structure of the subbands of double-occupied (left) and single-occupied 
(right) localized states for the electron occupation corresponding to metallic state in the 
"critical" region in the absence of magnetic field. 

Fig.2,b The structure of the subbands for high magnetic field driving the system to 
strongly localized regime. Note that the position of e m in the D° band is considered to be 
lifted with respect to the case H = (moreover, the form of the bands is also changed) 
while the position of fi in the D~ band is lowered. 

Fig. 3 Resistivity versus temperature for Si 5-doped GaAs-AlGaAs heterostructure of 
Ref. |13|] for B = 6 (full squares) and 8 Tesla (crosses) at n — 9.52 x 10 10 cm~ 2 . Solid lines 
- theoretical curves corresponding to our Eq.|TJ] calculated for e c — \i = 0.234meV(H = 
8T),0.219meV(H = 6T), e x = O.UmeV , g a 2 = 0.93meV-\ g = 0.12. 

Fig.4 Resistivity versus temperature for Si MOSFET of Ref. @ for B = 1.2 (full 
squares) and 1.4 (crosses) at n = 8.83 x 10 10 cm -2 . Solid lines - theoretical curves calcu- 
lated on the base of Eq.|l4| with the values of the parameters e c — \x = 0.147meV(H = 
l.4T),0.126meV(H = 1.2T). e 1 = O.lmeV, g a 2 = l^meV' 1 , g = 1.7. 
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